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ABSTRACT 

We investigate the differences in the M, - a relation derived recently by Ferrarese & Merritt and by Gebhardt 
and collaborators. The shallower slope found by the latter authors (3.75 vs. 4.8) is due partly to the use of a 
regression algorithm that ignores measurement errors, and partly to the value of the velocity dispersion adopted 
for a single galaxy, the Milky Way. A steeper relation is shown to provide a better fit to black hole masses derived 
from reverberation mapping studies. Combining the stellar dynamical, gas dynamical, and reverberation mapping 
mass estimates, we derive a best-fit relation M. = 1 .30(±0.36) x 10 8 M Q ((7 c /200 km s - 1 )4-72(±0.36) > 
Subject headings: black hole physics — galaxies: kinematics and dynamics — methods: data analysis 



1. INTRODUCTION 

Ferrarese & Merritt (2000; FM00; Paper I) demonstrated 
a tight correlation between the masses of supermassive black 
holes (BHs) and the velocity dispersions of their host bulges, 
M. oc a a , a = 4.8 ± 0.5. The scatter in the relation was found 
to be consistent with that expected on the basis of measure- 
ment errors alone; in other words, the underlying correlation 
between a and M. is essentially perfect. The relation is ap- 
parently so tight that it surpasses in predictive accuracy what 
can be achieved from detailed dynamical modelling of stel- 
lar kinematical data in most galaxies. As an example, FM00 
showed that the BH mass estimates of Magorrian et al. (1998), 
derived from ground-based optical observations, lie systemati- 
cally above the M. - a relation defined by galaxies with secure 
BH masses, some by as much as two orders of magnitude. 

The M. - a relation of Paper I was based on central velocity 
dispersions a c , corrected to an effective aperture of radius r e /8, 
with r e the half-light radius. Central velocity dispersions are 
easily measured and available for a large number of galaxies 
(Prugniel et al. 1997). An alternative form of the M m - a rela- 
tion was investigated by Gebhardt et al. (2000a; GOO) who used 
a e as the independent variable; a e was defined as the spatially- 
averaged, rms, line-of-sight stellar velocity within the effective 
radius r e . Computing <j e requires knowledge of the stellar ro- 
tation and velocity dispersion profiles at all radii within r e , as 
well as information about the inclination of the rotation axis 
with respect to the line of sight. These data are available for 
a smaller number of galaxies; on the other hand, <j e might be 
expected to reflect the depth of the stellar potential well more 
accurately than a c . 

The versions of the M, - a relation derived by FM00 and by 
GOO differ in two important ways. The latter authors found a 
significantly smaller slope (a = 3.75 ±0.3 vs. 4.8 ±0.5) as well 
as a greater vertical scatter - greater both in an absolute sense, 
and relative to measurement errors in M, . GOO estimated that 
approximately 40% of the scatter in M, about the mean line was 
intrinsic and the remainder due to measurement errors. FM00 
found no evidence for an intrinsic scatter in M.. 

The M. - a relation is currently our best guide to BH demo- 
graphics, and it is important to understand the source of these 

1 The error bars plotted in Fig. 2 of GOO do not always correspond to the values 
values. 
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differences. That is the goal of this paper. In addition to us- 
ing different measures of the velocity dispersion, FM00 and 
GOO analyzed different galaxy samples, and used different al- 
gorithms for fitting regression lines to the data. We find that 
regression algorithms that account correctly for errors in the 
measured variables always give a steeper slope than that found 
by GOO. We also show that the steeper relation derived by FM00 
provides a better fit to galaxies with BH masses computed by 
reverberation mapping. 



2. DATA 

Table 1 gives the data used here. The first 12 galaxies (Sam- 
ple 1) are "Sample A" from Paper I, consisting of those galaxies 
with published BH mass estimates that were deemed reliable - 
roughly speaking, galaxies in which the sphere of influence of 
the BH has been resolved. Five of these masses are derived 
from stellar kinematics and seven from gas dynamics. All of 
these galaxies were included in the GOO sample as well, with 
the exception of NGC 3115; for this galaxy, we assume a e = a c . 
The second part of Table 1 contains the additional 15 galaxies 
included by GOO (Sample 2). 1 Most of the BH mass estimates 
for these galaxies are based on unpublished STIS data. In addi- 
tion, GOO included M31 and NGC 1068, which were excluded 
from Paper I on the grounds that their BH masses were deemed 
unreliable. We computed distances for the GOO galaxies in Ta- 
ble 1 in the same way as in Paper I and corrected the BH masses 
accordingly. We also computed aperture-corrected central dis- 
persions cr c for the GOO galaxies. 

At this point, we are already in a position to test the idea, 
proposed by GOO, that the steeper slope of the M, - a c rela- 
tion in Paper I is due to spuriously high values of a c for the 
more nearby galaxies. This idea is rejected based on Figure 1, 
which shows that there is remarkably little difference on aver- 
age between a c and cr e . This is presumably due to the flatness 
of galaxy rotation and velocity dispersion profiles, and to the 
fact that even a c is measured on a large enough scale that it is 
essentially unaffected by the presence of the BH. The mean ra- 
tio of a e to er c is 1.01; the correlation coefficient of logcr^ vs 
log(T c is 0.97. 

listed in their Table 1 (e.g. NGC 4291 and NGC 5845). We used the tabulated 
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Table 1. BH Masses and Galaxy Velocity Dispersions 



G <il tixy 


Tvne 


Distance 


M. 










Sample 1 (FMOO) 






MW 


Sbl-II 


0.008±0.0009 


0.0295±0.0035 


100±20 


75 


11459 


E3 


30.3±4.0 


4.6±2.8 


312±41 


323 


N221 


cE2 


0.8±0.1 


0.039±0.009 


76±10 


75 


N3115 


so- 


9.8±0.6 


9.2±3.0 


278±36 




N3379 


El 


10.8±0.7 


1.35±0.73 


201±26 


206 


N4258 


SAB(s)bc 


7.2±0.3 


0.390±0.034 


138±18 


120 


N4261 


E2 


33.0±3.2 


5 4 +1 - 2 


290±38 


315 


N4342 


SO" 


16.7±1.0 


3.3+;-? 


261±34 


225 


N4374 


El 


18.7±1.2 


17+ 12 
1 '-6.7 


286±37 


296 


N4486 


EOpec 


16.7±1.Q 


35.7±10.2 


345±45 


375 


N6251 


E 


1Q4±1Q 


5.9±2.0 


297±39 


290 


N7052 


E 


66.1±6.4 


, ,+2.6 
J -'-1.5 


261±34 


266 






Sample 2 (GOO) 






N821 


E6 


24.7±2.5 


0.51±0.2 


196±26 


209 


N224 


Sb 


0.77±0.04 


0.35±0.25 


112±15 


160 


N1023 


SO 


10.7±0.8 


0.39+°-™ 


201±14 


205 


N1068 


Sb 


23.6±3.2 




149±19 


151 


N2778 


E 


23.3±3.4 


0.20^ 


171±22 


175 


N3377 


E5+ 


11.6±0.6 


03 +i - G 


131±17 


145 


N3384 


SB(s)0" 


11.9±0.9 


n 1 oc+0.072 
°' 185 -0.091 


151±20 


143 


N3608 


E2 


23.6±1.5 


1 13+ 1 - 44 




182 


N4291 


E 


26.9±4.1 


1 54+ 3 - 1 

1-0 -0.68 


269±35 


242 


N4473 


E5 


16.1±1.1 


I.O26+071 


188±25 


190 


N4564 


E 


14.9±1.2 


U - a '-0.17 


153±20 


162 


N4649 


E2 


17.3±1.3 




331±43 


375 


N4697 


E6 


11.9±0.8 


, 99 +o.io 


163±21 


177 


N5845 


E* 


28.5±4.2 


a - oz -0.72 


275±36 


234 


N7457 


SA(rs)0" 


13.5±1.4 




73±10 


67 



NOTE. - Type is revised Hubble type. Black hole masses are in 10 solar 
masses. Distances are in Mpc. Velocity dispersions are in km s -1 . 



However, we notice that a e and a c differ significantly for one 
particular galaxy, the Milky Way. GOO adopted a value of <j e = 
75 km s -1 for the Galaxy based on the velocity dispersion be- 
tween 50 and 500 arcsec (Kent 1992; Genzel et al. 2000). 
They apparently neglected to account for the contribution of 
the rotational velocity, which in the same region is 103±15 km 
s _1 (Kent 1992). More importantly, 500 arcsec corresponds to 
a projected radius of 20 pc at the Galactic center, more than 
two orders of magnitude smaller than the effective radius of 
the Galactic bulge (<~ 2.7 kpc; cf. Gilmore, King & van der 
Kruit 1990). The bulge velocity dispersion has been measured 
by several authors at various Galactocentric distances within 4 
kpc, all giving values between 75 and 110 km s" 1 , with a ten- 
dency for a to increase slowly toward the center (e.g. Kent 
1992; Tiede & Terndrup 1997, 1999; Minniti 1996; Cote 1999; 
Zhao et al. 1996). The rotational velocity in the inner 1.5 kpc 
is well approximated by a solid body curve with v ~ 65 - 87 
km s-'kpc -1 (e.g. Tiede & Terndrup 1997, 1999; Morrison & 
Harding 1993; Menzies 1990; Kinman, Feast & Lasker 1988). 
In view of these results, we question the choice of a e = 75 km 
s _1 for the Milky Way. FMOO adopted a c = 100 km s" 1 ; in what 
follows we will perform regression analyses assuming values 



of both 75 km s _1 and 100 km s-'for 



We will show that 



the slope of the M. - a e relation depends significantly on which 
value is used. 

GOO assumed constant errors on logM. and zero measure- 
ment errors in a e when carrying out their least-square fits. How- 
ever, ignoring measurement errors in the independent variable 
is well known to bias the slope downwards (e.g Jefferys 1980). 
Even when high signal-to-noise data are used, measurement er- 
rors on the velocity dispersions are easily at the 10% level (e.g. 




100 



200 



500 



ct (km/s) 



FIG. 1 . — Comparison of <j c (the central velocity dispersion) and a e (the rms 
velocity within one effective radius) for the galaxies in Table 1 . The solid line 
has a slope of one. 



van der Marel et al. 1994), and cannot be neglected. Unfortu- 
nately, the data used by GOO to compute a e are mostly unpub- 
lished, and the authors do not give error estimates in their paper. 
Therefore, in the regression analyses described below we will 
make various assumptions about the measurements uncertain- 
ties in a e . 

To understand how the different galaxy samples used by 
FMOO and by GOO may have affected their respective conclu- 
sions about the M.-a relation, we analyze Sample 1 (the 12 
galaxies from Paper I) and Sample 2 (the additional 15 galaxies 
from GOO) separately. The BH masses in Sample 2 are signifi- 
cantly less accurate than those in Sample 1, with an rms uncer- 
tainty of 0.28 dex, compared to 0.18 dex for Sample 1. We also 
present results from the analysis of the entire set of 27 galaxies, 
called the "combined sample" below. 

3. ANALYSIS 

We assume a relation of the form 



Yi = aXi + f3+ei 



(1) 



between the measured variables, where Y is logM. and X is 
either log a c or logcv Units are solar masses for M. and km 
s"'for a. The e, describe measurement errors as well as intrin- 
sic scatter in the relation, if any. A large number of regression 
algorithms are available for recovering estimates a and $ of the 
slope and intercept and their uncertainties, given (Xj, Yf) and 
their estimated errors (axi, cy;)- These algorithms differ in the 
degree of generality of the model that is assumed to underly the 
data. The following four algorithms were used here. 

Ordinary Least-Squares (OLS). All of the error is assumed 
to lie in the dependent variable (i.e. logM.) and the amplitude 
of the error is assumed to be the same from measurement to 
measurement. This is the algorithm adopted by GOO. We use 
the implementation G02CAF from the NAG subroutine library. 
The OLS estimator is biased if there are measurement errors 
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Table 2. Results of the Linear Regression Fits: logM. = alogcr c + f3 



Fit 




Combined Sample 






Samp] 


e 1 






Samp! 


e 2 






a 


a(a) 


P 


HP) 


a 


a(a) 


& 


HP) 


a 


Ha) 


P 


HP) 


OLS 


4.00 


0.35 


-1.10 


0.80 


4.39 


0.38 


1.94 


0.89 


3.28 


0.61 


0.47 


1.37 


GLS 


4.49 


0.13 


-2.18 


0.28 


4.71 


0.14 


2.64 


0.31 


3.32 


0.38 


0.40 


0.86 


ODR 


4.52 


0.36 


-2.31 


0.83 


4.46 


0.43 


2.06 


1.03 


4.20 


0.60 


1.66 


1.37 


BRS 


4.43 


0.39 


-2.08 


0.92 


4.81 


0.55 


2.92 


1.30 


3.75 


0.59 


0.59 


1.32 



Table 3. Results of the Linear Regression Fits: logM. = alogcTe + /3 



Fit 




a 


Combined Sample 

<r( Q ) a 


HP) 


a 


Sample 1 
a(a) fj 


HP) 


a 


Sample 2 

o(a) 


HP) 


OLS 




3.85 


0.30 


-0.75 


0.68 


3.86 


0.29 


-0.65 


0.67 


3.53 


0.55 


-0.13 


1.25 






(4.00)< 2 > (0.32) 


(-1.10) 


(0.74) 


(4.12) 


(0.38) 


(-1.30) 


(0.90) 










GLS 




3.90 


0.11 


-0.73 


0.23 


3.93 


0.12 


-0.74 


0.25 


3.58 


0.039 


-0.21 


0.089 






(4.47) 


(0.13) 


(-2.14) 


(0.28) 


(4.71) 


(0.14) 


(-2.64) 


(0.31) 










ODR 


5% 


3.90 


0.15 


-0.81 


0.35 


3.95 


0.16 


-0.83 


0.37 


3.96 


0.43 


-1.08 


0.99 






(4.33) 


(0.17) 


(-1.83) 


(0.39) 


(4.42) 


(0.19) 


(-1.96) 


(0.43) 












10% 


4.07 


0.24 


-1.20 


0.55 


3.98 


0.25 


-0.91 


0.59 


4.18 


0.51 


-1.58 


1.17 






(4.40) 


(0.27) 


(-2.01) 


(0.60) 


(4.42) 


(0.30) 


(-1.97) 


(0.70) 












15% 


4.17 


0.33 


-1.44 


0.75 


4.00 


0.36 


-0.96 


0.84 


4.27 


0.65 


-1.77 


1.49 






(4.45) 


(0.36) 


(-2.13) 


(0.83) 


(4.42) 


(0.43) 


(-2.00) 


(1.01) 












20% 


4.23 


0.43 


-1.59 


0.99 


4.02 


0.47 


-1.00 


1.11 


4.33 


0.82 


-1.92 


1.91 






(4.49) 


(0.47) 


(-2.23) 


(1.10) 


(4.44) 


(0.56) 


(-2.04) 


(1.35) 












25% 


4.27 


0.53 


-1.70 


1.24 


4.03 


0.58 


-1.04 


1.40 


4.40 


1.01 


-2.08 


2.40 






(4.53) 


(0.58) 


(-2.31) 


(1.37) 


(4.45) 


(0.71) 


(-2.06) 


(1.70) 










BRS 


5% 


3.90 


0.24 


-0.86 


0.55 


3.89 


0.22 


-0.73 


0.47 


3.60 


0.49 


-0.29 


1.10 






(4.05) 


(0.32) 


(-1.22) 


(0.74) 


(4.16) 


(0.41) 


(-1.40) 


(0.96) 












10% 


4.05 


0.27 


-1.20 


0.62 


4.00 


0.22 


-0.97 


0.47 


3.84 


0.57 


-0.83 


1.30 






(4.22) 


(0.36) 


(-1.61) 


(0.83) 


(4.29) 


(0.43) 


(-1.71) 


(1.01) 












15% 


4.33 


0.34 


-1.85 


0.79 


4.19 


0.24 


-1.41 


0.53 


4.32 


0.83 


-1.91 


1.91 






(4.54) 


(0.45) 


(-2.36) 


(1.04) 


(4.54) 


(0.48) 


(-2.28) 


(1.15) 












20% 


4.81 


0.50 


-2.95 


1.18 


4.49 


0.32 


-2.13 


0.74 


5.27 


1.65 


-4.04 


3.77 






(5.11) 


(0.65) 


(-3.64) 


(1.53) 


(4.94) 


(0.62) 


(-3.22) 


(1.48) 












25% 


5.64 


0.89 


-4.84 


2.08 


4.98 


0.53 


-3.25 


1.26 


7.42 


4.62 


-8.87 


10.5 






(6.11) 


(1.15) 


(-5.93) 


(2.70) 


(5.60) 


(0.93) 


(-4.75) 


(2.25) 











1 Percentages refer to assumed measurement errors in a e 

2 Values in parentheses used a e — 100 km s _1 for the Milky Way 



in the independent variable, or if the errors in the dependent 
variable vary from data point to data point (e.g. Jefferys 1980). 

General Least-Squares (GLS). All of the error is still as- 
sumed to reside in the dependent variable, but the amplitude 
of the error may vary from point to point. Press et al. (1989) 
implement this model in their routine fit, which we use here. 

Orthogonal Distance Regression (ODR). The underlying 
variables are assumed to lie exactly on a straight line, i.e. to 
have no intrinsic scatter, but the observed quantities are allowed 
to have measurement errors in both X and Y, which may differ 
from point to point. This model is incorporated in the routines 
f itexy of Press et al. (1989) and fv of Fasano & Vio (1988). 
We use the former routine here; the latter was found to give es- 
sentially identical results. The ODR estimator may be biased if 
the true variables exhibit intrinsic scatter about the linear rela- 
tion, in addition to measurement errors (e.g. Feigelson & Babu 
1992). 

Regression with Bivariate Errors and Intrinsic Scatter 
(BRS). As in ODR, the data are permitted to have measure- 
ment errors in both X and Y that differ from point to point. In 
addition, the underlying variables are allowed to have an in- 
trinsic scatter about the regression line. We use the routine 
BCES(F|X) of Akritas & Bershady (1996), the same routine 
used in Paper I. 



Tables 2 and 3 give estimates of the slope and intercept, a 
and $, and their uncertainties as computed by each of the four 
algorithms, using a c and a e as independent variables. Values in 
parentheses correspond to setting a e = 100 km s _1 for the Milky 
Way, as discussed above. The results are summarized below 
and in Figure 2. 

1 . Accounting for errors in one or both variables increases 
the slope of the relation, whether expressed in terms of a c or 
a e . Ignoring measurement errors biases the slope too low, for 
two reasons. The BH masses in Sample 2 are significantly more 
uncertain than those in Sample 1 and, as a group, exhibit a shal- 
lower slope (particularly when expressed in terms of a c ); rou- 
tines like OLS that weight all data equally therefore underes- 
timate the true slope. Second, ignoring measurement errors in 
the independent variable (log cr) always yields spuriously low 
slopes (e.g. Jefferys 1980). We find that the shallowest slope 
for every sample is returned by OLS, the routine used by GOO. 
All other algorithms give slopes in the range 4 < a < 5 for the 
combined sample. 

2. The slope inferred for the M, - <r e relation using ODR and 
BRS depends somewhat on the assumed errors in a e . Increas- 
ing the assumed error from 5% to 20% increases the BRS slope 
of the combined sample from 3.9 to 4.8. 
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Fig. 2. - Regression fits to the three galaxy samples, with a c or a e as independent variable. Sample 1: galaxies from Paper 
I. Sample 2: additional galaxies from GOO (also Table 1). Regression lines were computed with the BRS algorithm, assuming 
10% errors in a e - The fit to Sample 1 is repeated in the right two panels as a dashed line. 



3. Even when the appropriate fitting routines are used, the 
M, - a c relation tends to have a steeper slope than the M, - a e 
relation. This difference, however, is driven by one galaxy only: 
when a e for the Milky Way is increased from 75 km s" 1 (used 
by GOO) to a more appropriate - in our opinion - value of 100 
km s~\ both relations have a best-fit slope of <~ 4.5 ±0.5 for 
the combined sample (assuming a plausible 10% - 15% error 
on a e ). 

4. Adding the galaxies of Sample 2 (from GOO) has little im- 
pact on the results, as long as measurement errors are taken into 
account by the fitting routine and u e = a c = 100 km s" 1 is used 
for the Galaxy: the regression lines for Sample 1 (from Paper I) 
and for the combined sample are essentially the same. In other 
words, the BH masses added by GOO are too uncertain to sig- 
nificantly alter the fit determined by the galaxies from Sample 
1 alone. 

We conclude that the different slopes found by GOO and 
FM00 (3.75 vs 4.8) are due partly to the neglect of measure- 
ment errors by the former authors, and partly to the difference 
between a e and a c for a single data point, the Milky Way. If we 
use the more appropriate value of a e = 100 km s _1 for the Milky 
Way, and a plausible 10% - 15% error on a e , the M,-a c and 
M,-(T e relations have essentially the same slope, ~ 4.5. The 
data points added by GOO, based mostly on unpublished 



modelling of stellar kinematical data from STIS, appear to con- 
tain little information about the M, — a relation that was not 
already contained in the more accurate masses from Paper I. 

We next address the scatter in the M.-a relation. The \ 2 
merit function for a linear fit to data with errors in both vari- 
ables is 

2 i N (ji-aXi-P^ 2 

(e.g. Press et al. 1989), where, in our case, Y = logM. and 
X = log a. A good fit has x 2 < 1 . Since measurement uncer- 
tainties ax,i are not available for the a e , we computed x 2 only 
for the M, - a c relation. We are also interested in the absolute 
scatter in logM., which we define as 




We computed x 2 an d A» using the fits given by the BRS re- 
gression algorithm (Table 2). The results are: 
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FIG. 3. — Reverberation mapping masses for seven galaxies. Solid line is the 
M, -o c relation from Ferrarese & Merritt (2000); dashed line is the M. -a e 
relation from Gebhardt et al. (2000a). 



Sample 1 : x 2 = 0.74 A. = 0.26 
Sample 2: x 2 = 1.67 A. = 0.35 
Combined sample : x 2 = 1 -20 A. = 0.34 

The Sample 2 galaxies exhibit a larger scatter in logM. than the 
galaxies in Sample 1 (0.35 dex vs. 0.26 dex), consistent with 
their greater measurement uncertainties (Table 1). Furthermore 
the 12 galaxies from Paper I define a significantly tighter corre- 
lation, as measured by x 2 , than the 15 galaxies added by GOO 
(X 2 = 0.74 vs. x 2 = 1-67), or than the combined sample. Thus 
we confirm the conclusion of GOO that the scatter in their data 
about the best-fit linear relation exceeds that expected on the 
basis of measurement error alone. The large x 2 for Sample 2 
may indicate that the measurement uncertainties quoted by GOO 
are too small. 

4. REVERBERATION MAPPING MASSES 

A long-standing discrepancy exists between BH masses de- 
termined from stellar kinematics and from reverberation map- 
ping; the latter technique uses emission lines in active galactic 
nuclei (AGN) to probe the virial mass within the broad-line re- 
gion (Netzer & Peterson 1997). Since there are currently no 
galaxies with BH masses determined independently by the two 
techniques, any comparison must be statistical. The standard 
approach (e.g. Wandel 1999) has been to compare the average 
BH mass at a given bulge luminosity as computed from rever- 
beration mapping with the mass predicted by the Magorrian et 
al. (1998) relation; the latter is based on stellar kinematical 
data, mostly of low spatial resolution. The discrepancy is a fac- 
tor of ~ 20 in the sense that the reverberation-mapping masses 
are too low (Wandel 1999). This discrepancy has most often 
been attributed to some unspecified, systematic error in the re- 

2 The velocity dispersions plotted by Gebhardt et al. (2000b) are labelled a e even 



verberation mapping masses (e.g. Richstone et al. 1998; Faber 
1999; Ho 1999). 

FM00 showed that the Magorrian et al. masses fall systemat- 
ically above the M. - a c relation defined by galaxies with more 
secure BH mass estimates, some by as much as two orders of 
magnitude. The offset is strongly correlated with distance sug- 
gesting a systematic, resolution-dependent error in the Magor- 
rian et al. modelling. Much or all of the discrepancy with the 
reverberation mapping masses might therefore be due to sys- 
tematic errors in the Magorrian et al. masses, contrary to the 
usual assumption. Gebhardt et al. (2000b) tested this idea by 
plotting seven AGN BH masses against their M, -a e relation. 2 
We reproduce that plot here, as Figure 3. The fit is reasonable, 
although the points tend to scatter below the line. We also plot 
in Figure 3 the steeper M. - o c relation derived in Paper I (given 
here, in Table 2, as the BRS regression fit on a c for Sample 1). 
The steeper relation of Paper I is clearly a better fit. 

We stress that several of the AGN data points lie at the 
low-mass end of the distribution where the M. - a relation is 
strongly affected by uncertainties in the slope. Nevertheless, 
there would no longer appear to be any prima facie reason for 
believing that the reverberation mapping masses are systemati- 
cally in error. Furthermore the scatter in these masses about the 
M, - a relation appears to be comparable to that of the Sample 
2 data from GOO. We therefore carried out regression fits com- 
bining the reverberation mapping masses with Sample 1 and 
Sample 2, for a total of 34 galaxies. We assumed 50% mea- 
surement errors in the AGN M, and 15% errors in the a c . The 
results (using the BRS regression routine) were: 

a = 4.72 ±0.36, (3 = -2.75 ±0.82, (4) 

very close to the parameters derived in Paper I using Sample 1 
alone. This fit has x 2 = 1-H and A. = 0.35, about as good as 
obtained using the galaxies in Table 1 . 

5. SUMMARY 

We investigated the differences in the M, - a relation as de- 
rived by Ferrarese & Merritt (2000) and by Gebhardt et al. 
(2000a). The latter authors found a shallower slope (3.75 vs. 
4.8) and a greater vertical scatter, larger than expected on the 
basis of measurement errors alone. Three possible explanations 
for the differences were explored: different galaxy samples; 
different definitions of the velocity dispersion, central (<r c ) vs 
integrated (a e ); and different routines for carrying out the re- 
gression. The shallower slope of the GOO relation was found 
to be due partly to the use of a regression algorithm that does 
not account properly for measurement errors, and partly to the 
adoption of a value of a e for the Milky Way which is, in our 
opinion, implausibly low. The greater scatter seen by GOO is 
due to larger uncertainties associated with the additional BH 
masses included by them, mostly from unpublished STIS data. 
When measurement uncertainties are properly accounted for, 
the parameters of the best-fit relation derived from the com- 
bined samples of FM00 and GOO are essentially identical to 
those derived from the sample of FM00 alone. The steeper re- 
lation derived by FM00 also provides a better fit to BH masses 
obtained from reverberation mapping. A regression fit to the 
combined sample of 34 galaxies, including stellar dynamical, 
gas dynamical, and reverberation mapping masses, yields: 

though they are central values. 
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M.-o Relation 



M. = 1.30(±0.36)x 10 8 M Q . 

1 200kms _1 



4.72(±0.36) 



(5) 



The scientific implications of Equation (4) are discussed 
briefly by FMOO, and extensively in Merritt & Ferrarese (2000). 
This relation is essentially identical to the one derived in Pa- 
per I. We suggest that there is no longer any reason to assume, 
as a number of authors (Richstone et al. 1998; Faber 1999; 
Ho 1999) have done, that the reverberation mapping masses 
are less accurate than masses derived from stellar kinematics 
(Magorrian et al. 1998). 

We stress that the current sample of galaxies with reliable 
BH mass estimates is likely affected by severe selection biases, 
which are very difficult to quantify. Our results highlight the 



need for accurate BH masses if the M. - a relation is to be fur- 
ther refined. Only a handful of galaxies observed with STIS 
are likely to yield mass estimates as accurate as those already 
available for the galaxies in Paper I. Uncertainties in the re- 
verberation mapping masses are probably comparable to those 
obtained from HST data in most galaxies; however the number 
of galaxies with reverberation mapping masses is large (~ 35) 
and growing. Furthermore, many of these galaxies are in the 
critical, low mass range, 10 6 < M, < 10 8 M Q . An aggressive 
campaign to measure stellar velocity dispersions in AGN might 
be the best route toward refining the M, - a relation. 
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